Low mass dilepton radiation at RHIC 
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In this work we discuss the emission of low mass dilepton radiation from a hydrodynamic evolution 
model of Au-Au collisions and make comparisons with recent PHENIX measurements. The dilepton 
emission rates from the hadronic phase are treated at finite temperature and Baryon density and are 
completely constrained by broken chiral symmetry in a density expansion. The rates are expressed in 
terms of vacuum correlators which are measured in e + e~ annihilation, r decays and photo-reactions 
on nucleons and nuclei. We consider two possibilities for the hadronic phase: A chemical equilibrated 
an off equilibrium hadronic gas. We find that while chemical off-equilibrium helps explain part of 
the low mass (0.15 < M GeV < 0.7) enhancement seen in the data there is still a large discrepancy. 
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Dileptons and photons are a particularly interesting observable from heavy ion collisions since electromagnetic 
probes do not interact with the medium after they are produced and therefore carry information on the early stages of 
the evolution. In an ideal world one would hope to have an unambiguous signal coming from the quark gluon plasma 
in order to verify its existence and study its properties. In any collision, there is also a substantial contribution of 
dileptons coming from the hadronic phase. The hadronic emission dominates over the quark-gluon plasma signal and 
is therefore seen as a large background to one interested in studying the QGP. However, the hadronic phase does 
contain interesting physics in itself. An understanding of the resulting hadronic yields can provide crucial information 
on modifications to electromagnetic spectral functions due to chiral symmetry restoration. 

There is a long history of both experimental and theoretical work in electromagnetic probes, which we don't attempt 
to summarize here. Most recently there has been two experiments that have looked at dilepton emission in heavy- 
ion collisions. The recent NA60 experiment at the CERN SPS has measured the invariant mass spectrum and the 
transverse mass spectrum of low-mass dimuon pairs [l],[2|,[3| m I n ~I n collisions. It was seen that a large excess remained 
| after subtracting contributions from expected hadronic (the cocktail) decays. The remaining excess was examined 
by a number of groups 0, 0, S 0, HI, 0] and was interpreted as a combination of thermal partonic and hadronic 
(~\| ' contributions with modifications to the spectral function due to finite temperature and Baryon density. 

An experiment at RHIC, performed by the PHENIX collaboration PUHH, measured di-electron pairs in Au-Au 
ON ' collisions. They find a large excess above the cocktail in the mass region 0.15 < M (GeV) < 0.7 for more central 
collisions (N part > 250). For most central collisions (N part ~ 325) the yield in this mass region is increased by a 
' factor of ~ 8 above the cocktail. 

Quickly looking at the two most recent data sets from NA60 and PHENIX one might conclude that the data are 
inconsistent with each other. This is not necessarily the case as the two detectors have very complicated acceptances 
which distort the resulting yields. Also the space-time evolution and the resulting chemistry is different in In-In and 
Au-Au collisions. Therefore, before any conclusions can be drawn, it is necessary to calculate the resulting yields one 
expects from standard rate equations at RHIC. 

In this work we calculate the resulting di-electron yields from Au-Au collisions at RHIC energies. In the next section 
we summarize the rate equations used in both the partonic and hadronic phases. Then we describe the hydrodynamic 
evolution model tuned to RHIC collisions. We consider two scenarios in this work: A chemically equilibrated hadronic 
phase as well as a chemical off-equilibrium hadronic phase. The resulting yields are then compared to the recent 
PHENIX data. 



II. DILEPTON RATES 



In this section we summarize the dilepton rates used in the analysis. We should note that these are the same rates 
used in a previous analysis of the NA60 data 4]. For the partonic contribution above a critical temperature T c « 180 
MeV we use the standard leading order qq result for massless quarks: 



q=u,d,s 

where Nc is the number of colors, e q the charge of the quarks, and n± = l/(e^ 0± l^^/ 2T + 1). 
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Below T c we use the rate equations presented in [12|, |l3|, for a hadronic gas at finite temperature and Baryon 
density which are constrained entirely by broken chiral symmetry: 

f q = ^$TT^7? (1 + 2 -f ){l ~ ^) 1/2 h3, 2 Imn y (, 2 ) +±f d-Wf (*,*) + J dNW-( q , P )] (2) 

where da and dN are the appropriate phase space factors for Mesons and nucleons respectively as outlined in 0]. 
The term in square brackets is obtained after keeping terms to first order in an expansion of the thermal structure 
function in both Meson and nucleon density. Through the use of three-flavor chiral reduction formulas the terms 
~Wf (q, k) and Wjy(g,p) can be expressed in terms of vector and axial spectral densities which are measured in e + e~ 
annihilation, t decays and photo-reactions on nucleons and nuclei. 

In figure [1] we show the dilepton rates at two points selected so that one is representative of the beginning of 
the hadronic evolution and another towards the end. The left figure shows the rates for T=180 MeV and //b=27 
MeV which is the conditions at chemical freeezout. The curve labeled QGP is the leading order contribution from 
qq annihilation. The term labeled 'Vac' corresponds to contribution to zeroth order in the density expansion. The 
contribution from the term to first order in the density expansion in Mesons is separated into two parts: those from 
the vector spectral function (Hv) and those from the axial (If^). We should note for clarity that even though we 
show the QGP contribution in Fig. Q] at a temperature of 90 MeV it only exists in our analysis above the critical 
temperature. 

The main contribution comes from Ha which gives a strong enhancement in the rates below the two-pion threshold. 
The mechanism for this can be qualitatively thought of as X — > 7r + e + +e~. It is only qualitative because the physical 
kinetic processes are mixed by the virial expansion [l2j ■ 

The last piece 'Nucl' is the contribution from finite Baryon densities. The chemistry at RHIC does not initially 
support a Baryon rich environment in contrast with fixed target experiments. Therefore, at chemical freezeout the 
Baryon density is low and will not significantly modify the rates. This situation is shown in the left of fig. [T] Even 
though the net Baryon number is small at RHIC, the total Baryon number is not. This observation is taken into 
account by imposing conservation of both Baryon and anti-Baryon number in the hadro-chemical evolution and will 
be discussed at length in the following section. At lower temperatures we find a larger Baryon chemical potential but 
due to the Boltzmann suppression the rates are still largely unaffected by Nucleons as seen in the right of fig. [T] 
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FIG. 1: Dilepton rates used in this work. The QGP is the leading order qq annihilation result. The curve labeled 'Vac' is the 
first term in eq. [2] and corresponds to the resonance gas contribution. The curves labeled II a and II v correspond respectively to 
the axial- vector and vector contributions in the density expansion of eq.[2] The curve labeled 'Nucl' is the nucleon contribution. 
The bold solid curve is the total hadronic rate. Left: T=180 MeV fx B =27 MeV Right: T=90 MeV fi B = 428 MeV, ^=90 
MeV, ^=220 MeV. 



III. EVOLUTION 



In order to see if the rates given in the previous section reproduce the experimentally observed yield, the rates 
must be convoluted over the full space-time history of the collision region having widely varying temperature, Baryon 
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Parameter 


Value 




0.05c 


C QGP 


0.33c 


Tc 


180 MeV 


T/.o. 


120 MeV (90 MeV) 


TO 


0.6 fm/c 


n B /s 


0.004 


C s 


15 




0.06 


b 


3.0 fm 



TABLE I: Parameters used in the hydrodynamic simulation of Au-Au collisions at y/s = 200 GeV. A freezeout temperature of 
120 MeV is used for the hadronic phase in chemical equilibrium while T = 90 MeV is used for off-equilibrium. 

density, chemical potential and flow velocity. The evolution of the collision region is treated with a hydrodynamical 
model based on reference [i~o| . The model is a 2+1 dimensional boost invariant simulation with an ideal gas equation 
of state (p = |e) in the QGP phase and an ideal gas of Meson and Baryon resonances in the hadronic phase. 

The initial condition is set by the entropy in the transverse plane according to the distribution of participants for 
Au-Au collisions. One parameter, (C s =entropy per unit spatial rapidity), is adjusted to set the initial temperature 
and total particle yield. A second parameter (C„ B =Baryon number per unit spatial rapidity), is adjusted to fit the 
net yield of protons. The parameters used in the hydrodynamic simulation is summarized in table HI 

For the hadronic phase we consider two different scenarios. A chemical equilibrated phase (i.e., vanishing chemical 
potential for all particle species except Baryons) or a hadronic phase out of chemical equilibrium. In both cases the 
relationship between pressure and energy density remain about the same (see Appendix A) and therefore the hydro- 
dynamic solution remains the same. However, the temperature as a function of energy density changes significantly. 
The hadronic phase cools much quicker out of equilibrium. 

As discussed in detail in the Appendix thermal equilibrium is maintained for a much longer time period in the 
evolution compared to chemical equilibrium in the hadronic phase [l6| . This is due to the elastic scattering cross 
section being much larger then the inelastic cross section. Elastic collisions do not change the net quantum numbers 
of particles. For example the main hadronic reactions, e.g., tttt — » p — > t:t:,ttK — > K* — > ■nK,TiN — > A — > ttN do 
not change the net yield of Pions, Kaons and nucleons. The inelastic collisions which do change the yield of hadronic 
species take place on time scales longer then the collision time. In order to take this into account in hydrodynamic 
simulations which implicitly assume thermal and chemical equilibrium, effective chemical potentials are introduced 
which keep the particle yields constant throughout the evolution. Our analysis will therefore consist of two different 
scenarios. A chemical equilibrated phase where the chemical potential /i vanishes for all species except Baryons and 
a chemical off-equilibrium phase where rij/s is conserved throughout the evolution, ni corresponds to the number 
density of any hadronic specie listed in the appendix and s is the entropy density. 

There is one further point worth making concerning the role of conserved species. Rapp [i"7j showed the importance 
played by conservation of anti-Baryon number on the abundances of other species. Without anti-Baryon number 
conservation the Meson chemical potentials tend to remain small. However, when separately imposing anti-Baryon 
conservation the Meson chemical potentials become quite large. In our work fi n w 90 MeV and [Ik ~ 200 MeV at 
T » 100 MeV. The reason for this is the large amount of entropy stored in BB pairs. This results in a smaller amount 
of entropy per Pion therefore requiring a larger /i^ in order to keep n n /s constant. 

In figure [2] the chemical potential for Pions, Kaons, Nucleons and anti-Nucleons is shown as a function of temper- 
ature. Details of the calculation are discussed in the appendix. We note that in the chemically equilibrated scenario 
there is only one chemical potential, /is, corresponding to conservation of Baryon number. 

In figure [3] we show the evolution of temperature at the center of the collision region as a function of proper time. 
For the hadronic phase we show two trajectories, one for the chemical equilibrated and the other for the off equilibrium 
hadronic gas. We find that the hadronic phase cools quicker in the off-equilibrium (finite fi) scenario. 

In figure 0] the hydrodynamic solution for Semi-Central collisions is shown. The three contours correspond to 
the phase transition from the QGP to mixed phase, from the mixed to hadronic phase and the last contour is the 
space-time location of kinetic freezeout. 

IV. RESULTS & DISCUSSION 

We now show the results after evolving the rate equations discussed above through the hydrodynamic evolution. 
We show the momentum integrated rates dN/dM for a slice in rapidity \y\ < 0.7 corresponding to the PHENIX 
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FIG. 2: Evolution of chemical potentials as a function of temperature. 
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FIG. 3: Temperature at the center of the hydrodynamic simulation of Au-Au collisions at RHIC as a function of r at an impact 
parameter of b=3 fm. The curve labeled ('with denotes the trajectory when chemical off-equilibrium is used. Freezeout 
occurs in the center at r ~ 14 fm/c. 

detector acceptance in figure 

The partonic contribution does not become dominant until after the <f> mass. In the low mass region of interest 
here it is sub-dominant by more than a factor of two. The overall yields from the two different scenarios used in the 
hadronic phase are comparable. This should be expected. Even though the fireball temperature drops quicker in the 
off-equilibrium scenario the dilepton rates are enhanced by additional fugacity factors. Near the p peak at M=770 
MeV this corresponds to a factor of z 2 where z = e M,r / T which compensates from the reduction due to the Boltzmann 
factor. Below the p peak the yields are enhanced by a factor of sa z 3 . The extra factor of z comes from the Meson 
phase space factor in eqn[21 

This difference in fugacity factors near and below the p mass region is what modifies the shape of the hadronic 
yields when including off-equilibrium effects. The larger Pion and Kaon densities in the off-equilibrium hadronic phase 
shuffles strength from the rho peak into the low mass region. At M « 0.4 MeV we see about a 50% enhancement in 
the yields. 

In figure [5] we show the yields from the chemical off-equilibirum scenario after being evolved through the full 
detector cuts at PHENIX. The overall yield is determined by rescaling the results by the number of participants 
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FIG. 4: The hydro dynamic solution for Semi-Central (b=3 fm) Au-Au collisions at RHIC assuming chemical equilibrium. The 
solid lines show contours of constant energy density as a function of y or x at x=0 or y=0 respectively. The contour values are 
for e q — 1.8 GeV/fm 3 , = 0.65 GeV/fm 3 and for ef, a . = 0.07 GeV which is corresponds to a kinetic freezeout temperature of 
T = 130 MeV or T = 90 MeV in and out of equilibrium respectively. The two curves correspond to slices in the x and y planes 
of the simulation. 
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FIG. 5: Total integrated yield for RHIC showing the partonic contribution and the hadronic contribution with and without 
chemical freezeout. 



corresponding to the minimum Bias data set. In going from our central collision to min. bias this corresponds to a 
factor of ~ N par t(min.bias) /Npart(central) = 109/325. Ideally we should compare our yields to the 10% centrality 
class, which will be done when the data is published. Even though chemical off-equilibrium helped explain part of the 
excess di-electrons in the low mass region there is still a large part of the yield that is unexplained. In an attempt to 
reconcile this we have tried one further scenario. A chemically super-saturated phase throughout the entire hadronic 
lifetime. A similar analysis was performed in the framework of a Boltzmann type transport model by [20| . In this 
model it is assumed that the hadronic phase starts as a supersaturated pion gas with large effective potentials. 

We perform this calculation schematically in order to get a qualitative picture of the result. We use the same 
evolution model as was used for the chemically equilibrated phase but introduce a pion chemical potential of \i v = 50 
MeV in the rates which remains constant during the entire phase. In a more realistic picture one would also need 
to include this chemical potential in the hydrodynamic equation of state as well. This would cause the hadronic 
phase to cool much quicker. Since we are not doing this we obviously expect to overestimate the yields. Yet it is still 
interesting to see how the inclusion of such a chemical potential changes the shape of the spectrum. This result is 
shown in fig. [lV]by the curve labeled 'SS'. 

We should mention that R. Rapp [H GJ] has also performed a similar analysis which included thermal partonic 
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FIG. 6: Total integrated yield for RHIC showing the sum of partonic contributions and the hadronic contribution with chemical 
freezeout after a rescaling by the number of participants to fix the overall yield. 
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FIG. 7: Total integrated yield for RHIC showing the partonic contribution and the hadronic contribution with and without 
chemical freezeout. Also shown is the yield using a schematic supersaturated pion gas. 



and hadronic contributions. The hadronic contribution contains medium modifications in the form of self-energy 
corrections to the p, oj and (f> propagators from resonant interactions with the surrounding Mesons and Baryons as 
well as corrections from polarization from the pion cloud. Even though the form of the medium modifications differ 
in the two calculations the resulting yields are comparable. 



V. CONCLUSIONS 



In conclusion we have calculated the yields of di-electron pairs expected from thermal sources from Au-Au 
collisions at RHIC. We have considered a chemically equilibrated and a chemically off-equilibrium hadronic 
phase. Even though off-equilibrium effects help explain part of the low mass enhancement seen in the data 
there is still a large discrepancy. We have also shown the results from a schematic super-saturated pion gas, which 
when the correct normalization is taken into account, will probably have only a minimal effect in explaining the excess. 



7 



Acknowledgments. 

We would like to thank Alberica Toia for running our Yields through the PHENIX acceptance. KD would like to 
thank D. Teaney for help implementing chemical freezeout. This work was partially supported by the US-DOE 
grants DE-FG02-88ER40388 and DE-FG03-97ER4014. 

APPENDIX A: EQUATION OF STATE 

In this appendix we discuss some of the details in generating the equation of state used in the above analysis. The 
work below is based on previous calculations for SPS done by D. Teaney The necessity to introduce a separate 
chemical freezeout temperature arises because of the two different time scales present in heavy ion collisions: r t h and 
t c /j. The first timescale is the thermal equilibration time, r t h ~ 2 fm/c at T « 160 MeV and the second is the chemical 
equilibration time t c h ~ 200 fm/c. The time scale of a typical heavy-ion collision is on the order of Thic ~ 10 fm/c. 
In using a hydro dynamic description it is implicitly assumed that Tth *C tmc *C T c h- 

One therefore has the following picture: At the critical temperature hadronization occurs and the chemical com- 
position of all species is fixed. Since the time scale for chemical equilibrium is much longer then the lifetime of 
the collision the density of the particle species remain fixed (causing the particles to develop chemical potentials as 
the temperature decreases). The hadronic fluid remains in kinetic equilibrium due to elastic collisions until a final 
freezeout temperature when the hydrodynamic simulation terminates and we assume free streaming particles remain 
in the final state. 

The equation of state during the hadronic phase is constructed by assuming additional conservation laws for the 
following currents: 

B, 8,1,8 + s, 7T, K, 77, u, Y, Y, S, S, rf, fl, ft, 0, B (Al) 

We note that the conservation of anti-Baryons was not considered in (l6| but was taken into account due to its 
importance as pointed out in [13]. 

The procedure used to generate the equation of state is as follows. First at T c the values of fii and fi s are fixed 
so that s/tib — 250, nj = and n s = 0. Then calculate nu/s for the remaining 14 densities using their chemical 
equilibrated (/i = 0) values. Then in small increments the temperature is lowered and the chemical potentials adjusted 
in order to keep n^/s constant. 

We now discuss the results of this procedure. In figure [5] we show the pressure and speed of sound as a function of 
energy density. We see that chemical equilibrium does not change the result by much. We can therefore assume that 
the hydrodynamic solution does not change and use the same result for both the equilibrated and un-equilibrated 
phase making sure to use the appropriate temperature when calculating the rates and freezeout. 

The temperature as a function of energy density is shown in figure [51 We see that by introducing chemical freezeout 
the simulation cools much quicker. This will cause freezeout to occur earlier in the simulation and lead to shorter 
lifetimes. The smaller space-time volume will be compensated for in the yields by the fugacities in the rates. We also 
show the entropy density as a function of temperature in order to check our numerics. This should be the same both 
in and out of equilibrium, which is indeed the case as shown in figure O 

The final result for the chemical potentials as a function of temperature is shown in figure [2] A discussion of the 
results is kept to the text as it is relevant to the results. 
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FIG. 8: The pressure (left) and sound speed squared (right) as a function of energy density with and without chemical freezeout 
at RHIC (s/riB = 250). 
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FIG. 9: Entropy density as a function of energy density (left) and energy density as a function of temperature (right) with and 
without chemical freezeout at RHIC (s/ub = 250). 
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